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Abstract 

Retrievals of aerosol optical depth (AOD) from the Cloud-Aerosol Lidar with Orthogo- 
nal Polarization (CALIOP) satellite sensor require the assumption of an extinction-to- 
backscatter ratio, also known as the lidar ratio. This paper evaluates a new method 
5 to calculate the lidar ratio of sea spray aerosol using two independent sources: the 
AOD from the Synergized Optical Depth of Aerosols (SODA) algorithm and the in- 
tegrated attenuated backscatter from CALIOP With this method, the particulate lidar 
ratio can be derived for individual CALIOP retrievals in single aerosol layer columns 
over the ocean. Global analyses are carried out using CALIOP level 2, 5 km sea spray 
io aerosol layer products and the collocated SODA nighttime data from December 2007 
to December 2009. The global mean lidar ratio for sea spray aerosols was found to 
be 26 sr, roughly 30% higher than the current value prescribed by CALIOP standard 
retrieval algorithm. Data analysis also showed considerable spatiotemporal variability 
in the calculated lidar ratio over the remote oceans. The calculated aerosol lidar ratios 
is are shown to be inversely related to the mean ocean surface wind speed: increase in 
ocean surface wind speed (L/ 10 ) from 0 to > 15 ms -1 reduces the mean lidar ratios for 
sea spray particles from 32 sr (for 0 < L / 10 < 4 ms -1 ) to 22 sr (for U :Q > 15 ms -1 ). Such 
changes in the lidar ratio are expected to have a corresponding effect on the sea spray 
AOD. The outcomes of this study are relevant for future improvements of the SODA and 
20 CALIOP operational product and could lead to more accurate retrievals of sea spray 
AOD. 

1 Introduction 

Radiative forcing by sea spray aerosol (SSA) comprises a significant portion of the 
global energy budget. Studies have shown that SSA contributes an aerosol optical 
25 depth (AOD) of approximately 0.15 and cloud condensation nuclei (CCN) of about 
60 cm -3 (Kaufman et al., 2002; Lewis and Schwartz, 2004). Thus, SSA is an important 
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natural aerosol affecting both direct (i.e., extinction of solar radiation via scattering and 
absorption) and indirect (i.e., cloud lifetime and frequency) radiative forcing of climate. 
As SSA contributes substantially to the preindustrial, natural background and provides 
the base line on top of which anthropogenic forcing should be quantified, it is very 
5 important to properly characterise SSA burden and its spatiotemporal distribution. The 
representation of SSA in climate models was shown to strongly influence the predicted 
impact of anthropogenic aerosols on climate via direct and indirect effects (Ghan et al., 
2001 ; Hoose et al., 2009; Meskhidze et al., 201 1 ; Wang and Penner, 2009; Westervelt 
et al., 2012). Furthermore, uncertainties in preindustrial aerosol, of which sea spray 
io is a part of, were shown to contribute up to 45% in the overall variance in radiative 
forcing (Carslaw et al., 2013) and contribute the largest variance in model estimates of 
the aerosol indirect effect (Ghan et al., 2013). 

Aerosols over the remote oceans consist mainly of a mixture of marine (e.g., sea 
spray particles and sulfates from the oxidation of dimethylsulfide), natural continen- 
15 tal (e.g. mineral dust and biomass burning), and human-induced pollution (Andreae, 
2007). Therefore, knowing horizontal and vertical distribution, as well as speciation 
of aerosols becomes extremely important for the correct quantification of SSA radia- 
tive properties. The last decade has produced a large body of information regarding 
the sources and composition of marine aerosols, resulting in a reassessment of the 
20 complex role that sea spray particles play in climate and various geophysical phe- 
nomena. Passive satellite instruments like the Sea-Viewing Wide Field-of-view Sensor 
(SeaWiFS), the MODerate resolution Imaging Spectroradiometer (MODIS), and the 
Multi-angle Imaging Spectroradiometer (MISR), as well as the ground-based AErosol 
RObotic NETwork (AERONET) have contributed immensely to quantitative character- 
25 istics of SSA in terms of AOD (the column integrated aerosol extinction), size distri- 
bution information, and spectral optical properties. Although passive instruments have 
been useful for developing a basic picture of SSA distribution, they supply limited in- 
formation on aerosol speciation and very little data related to aerosol distribution in 
the vertical column. The introduction of the Cloud-Aerosol Lidar with Orthogonal Po- 
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larization (CALIOP) onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite 
Observations (CALIPSO) platform has eliminated some of the assumptions made by 
the passive instruments and has provided a more complete picture of global aerosol 
distribution desirable by climate scientists. However, CALIOP is an elastic backscatter 
5 lidar with no molecular filtering capability and therefore requires the assumption of an 
extinction-to-backscatter ratio, also known as the lidar ratio, to infer extinction from at- 
tenuated backscatter measurements. Depending on the microphysical properties of the 
aerosol, the lidar ratio can have a wide range of values and therefore a straightforward 
a priori solution within some reasonable uncertainty range is generally unobtainable 
10 without various assumptions or constraints. Theoretical calculations for the lidar ratio 
can be performed, if the physicochemical properties and the size distribution of the par- 
ticles at the different heights in the vertical column are known; although, the fulfillment 
of these requirements will make the lidar measurements unnecessary (Ackermann, 
1998). The typical solution to this problem is to assign a vertically independent lidar ra- 
15 tio to aerosol retrievals that fit a specific aerosol model as outlined in Omar et al. (2009). 

Since the uncertainty in the lidar ratio can significantly affect the accuracy of the 
aerosol extinction retrieval (see a detailed discussion below), lidar ratios have been 
constrained by numerous approaches. However, SSA size distribution, chemical com- 
position and refractive index can change significantly with ocean surface wind speed 
20 (L/ 10 ), relative humidity (RH), temperature, salinity and chemical/biological composition 
of surface sea water (de Leeuw et al., 2011; Lewis and Schwartz, 2004). Because of 
this, large disagreement exists in the literature regarding the value of maritime aerosol 
lidar ratio (S p ). For example, lidar measurements of Ansmann et al. (2001) over the 
North Atlantic showed S p = 24±5sr whereas measurements using a nighttime lidar 
25 measurements at a horizontal orientation off the northern coast of Queensland, Aus- 
tralia showed maritime aerosol lidar ratios as high as S p = 39 ± 5 (Young et al., 1993). 
Using the data from AERONET oceanic sites (Cattrall et al., 2005) derived a lidar ratio 
of 28 ± 5 sr, the value that compared well with literature averaged value of S p = 29 ± 5 sr 
(for 490 < A < 550 nm) for maritime aerosols. Passive techniques have also been used 
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to derive the lidar ratio using an alternative definition of S p as a function of single scat- 
tering albedo and the scattering phase function near 180° (Breon, 2013). Using the 
multi-directional measurements of the solar radiation from the polarization sensitive 
passive radiometer POLDER, typical values for clean marine aerosol S p were derived 
5 to be 25 sr at 532 nm (Breon, 201 3). The lidar ratio of 20 ±6 sr (at 532 nm) was selected 
for the CALIOP retrieval algorithm based on parameters measured during the Shore- 
line Environmental Aerosol Study (SEAS) experiment (Masonis et al., 2003; Omar 
et al., 2009). The SEAS experiment reports a particulate lidar ratio S p = 25.4 ± 3.5 sr at 
532 nm based on the optical size measurements of sea spray aerosol and a modeled 
io value of S p = 20.3 sr (Masonis et al., 2003). Although the S p value used in the CALIOP 
marine aerosol model is the same as the one derived using an average SSA size distri- 
bution measured on the beach (downwind of an offshore reef), modeling studies show 
a wide range of S p values (from 10 to 90 sr) depending on particle size (Masonis et al., 
2003). Therefore, as size distribution (and chemical composition) of SSA may vary over 
is the oceans, a constant lidar ratio used in CALIOP algorithms may lead to erroneous 
retrievals of AOD. 

In this study, we present a new method for deriving lidar ratios for individual CALIOP 
retrievals of single aerosol layer columns over the ocean. To estimate S p for a strictly de- 
fined subset of CALIPSO data we have used the Synergized Optical Depth of Aerosols 
20 (SODA) product (described in Sect. 2.2). The S p values are calculated as a correction to 
achieve the best agreement between SODA and CALIPSO sea spray AOD values. Us- 
ing CALIPSO level 2 aerosol layer data for years 2007 to 2010, we have created a 3yr 
averaged climatology of clean marine aerosol lidar ratio over the globe. Analyses were 
also carried out to assess dependence of S p values on wind speed and estimate pos- 
25 sible error sources in our calculations. It should be noted that other lidars, such as the 
High Spectral Resolution Lidar (HSRL), are capable of measuring aerosol backscatter 
and extinction parameters independently and therefore do not require the lidar ratio to 
be prescribed (Piironen and Eloranta, 1994). Despite this, currently CALIPSO is the 
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only lidar that provides aerosol data at the vast spatiotemporal resolution required for 
global climate model comparison. 


2 Methods 

2.1 CALIPSO satellite 

The CALIPSO mission was launched on April 28, 2006. CALIPSO has been able to 
provide the scientific community with vertically resolved measurements of both aerosol 
and cloud optical properties like depolarization ratio (a measure of particle spheric- 
ity), AOD, and ice/water phase since June 2006. The CALIPSO payload includes 
a high-powered digital camera, an infrared radiometer, and a two-wavelength (532 and 
1064nm), near nadir, polarization sensitive, elastic backscatter lidar, CALIOP. 

The level 1 data algorithms are responsible for the geolocation and range deter- 
mination of the satellite and produce profiles of attenuated backscatter coefficients. 
Data in this work were obtained from the 5 km, level 2 operational products version 
3.01. Level 2 products have undergone various processing algorithms from the Se- 
lective Iterated BoundarY Locator (SIBYL), the Scene Classification Algorithm (SCA), 
and the Hybrid Extinction Retrieval Algorithm (HERA). First, SIBYL identifies layers, 
then the SCA identifies the type of feature (i.e., aerosol or cloud) and the subtype 
(i.e., aerosol type, ice/water phase), and finally, the HERA generates extinction pro- 
files for the feature. The theoretical basis of the algorithm can be found online at 
www-calipso.larc.nasa.gov/resources/project_documentation.php. 

2.2 Synergized Optical Depth of Aerosols (SODA) 

CloudSat was launched in 2006 with CALIPSO and was positioned in sun-synchronous 
orbit as part of the A-Train satellite constellation. CloudSat and CALIPSO have paved 
the way for new multi-sensor data products like SODA to be developed. The main 
instrument on CloudSat is the Cloud Profiling Radar (CPR), a nearly nadir looking 
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(0.16°) 94 GHz (« 3 mm; W-band) radar. The CPR, like CALIOP, can retrieve informa- 
tion on hydrometeor microphysical properties at different heights in a vertical column. 
The CPR signal is mostly attenuated by water vapor; however, for cloud free regions 
over the ocean, the CPR data can be used to retrieve AOD. A method developed by 
s Josset et al. (2008) and later expanded by Josset et al. (2010a) uses a combination of 
CALIOP and CPR measurements of the ocean surface reflectance to derive AOD. The 
design of SODA utilises the ratio of the ocean surface scattering cross section of the 
radar over the lidar to infer column optical depth for non-cloudy atmospheric columns. 
The radar signal attenuates mostly due to water vapor, while the lidar signal weakens 
io mostly due to aerosols. Therefore, once the radar signal is corrected for attenuation by 
water vapor and oxygen, the change in the radar-to-lidar signal ratio is directly related 
to aerosol abundance (Josset et al., 2008, 2010a). Therefore, by using observations 
from two different sensors, SODA can eliminate uncertainties induced by the CALIOP 
aerosol extinction algorithm over oceans. SODA AODs have been shown to be in very 
is good agreement with MODIS AOD retrievals (Josset et al., 2008). More detailed de- 
scription of the SODA technique and its application is given in Josset et al. (2008, 
2010a, b, 2011, 2012) 

2.3 Lidar ratio definition 

One of the biggest advantages of the SODA product is that it removes the depen- 
20 dence of the prescribed lidar ratio while still utilizing the active sensors to retrieve 
an AOD, thereby providing a means for independent evaluation of the lidar ratio. In 
the current study we use Eq. (4) from Josset et al. (2011) to estimate lidar ratio from 
CloudSat/CALlOP measurements of AOD values. Following Fernald et al. (1972), the 
particulate two-way transmittance at height Z can be written as: 

25 f 2 (Z) = 0- 2S p/o^pW dz (1) 
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where the lidar ratio at height Z can be defined as the ratio of the particulate extinction 
to backscatter (S p = Differentiating Eq. (1) with respect to vertical coordinate (z) 
gives the particulate backscatter at height Z: 


5 


P P (Z) = - 


i d r 2 (Z) 

2S p 7" 2 (Z) dZ 


( 2 ) 


Since atmospheric constituents (molecules and different particle types) can inter- 
act with the lidar beam at different heights, the lidar ratio using remotely sensed data 
cannot be uniquely defined for a given atmospheric column. However, the lidar ratio is 
a particle intensive property (i.e., dependent on particle type and not on the amount). 
io So, if we assume that there is only a single type of aerosol and it is homogeneously 
distributed throughout the atmospheric column and that molecular scattering is suffi- 
ciently removed by the CALIOP level 2 algorithms, then the column lidar ratio (S p ) can 
be expressed as the ratio of the particulate column integrated extinction (r^ = AOD) to 

the attenuated backscatter (T p ). With these assumptions in mind, integration of Eq. (2) 
is with respect to vertical coordinate gives the particulate lidar ratio as: 




/^/? p (z)r p 2 (z)dz 


(3) 


In Eq. (3) if we substitute the definition for two-way transmittance as 7" p = e~ 2Tp , 
the total particulate attenuated backscatter signal retrieved by the lidar as r p = 

20 /q /? p (z)7" 2 (z)dz, and consider that 7" p (0) = 1 , we get an equation for a columnar partic- 
ulate lidar ratio as: 



1 - e 2lp 


(4) 
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This equation allows us to calculate SSA lidar ratio from two independent sources: 
the AOD (i.e., Tp) from SODA and the integrated attenuated backscatter (l~ p ) from 

CALIOP. It should be noted that CALIOP estimation of r p is difficult for layers that 
are not bounded by clear air (Vaughan et al., 2004) and therefore require carefully de- 
signed data screening algorithms. In Sect. 4 we carry out an error analysis to verify 
that uncertainties in r p have a minimal effect on the retrieved lidar ratio. 

2.4 Data selection method 

As different aerosol sub-types have different lidar ratios, application of Eq. (4) to 
episodes when aerosols other than sea spray are present in the atmospheric column 
may lead to erroneous results for the calculated S p . To minimise the contamination of 

sea spray AOD and therefore S p by aerosol types other than SSA (e.g., anthropogenic 
pollution, biomass burning, and dust), we developed a strict scene selection algorithm. 
The algorithm first uses the feature classification flags in the CALIOP aerosol layer 
product. We start with clean marine aerosol that is identified based on surface type 
(as determined by the location of the satellite), total integrated attenuated backscat- 
ter y' > 0.01 km -1 sr -1 , and volume depolarization ratio S' < 0.05 (Omar et al., 2009). 
As multiple types of aerosols can be found within retrieved vertical profiles (e.g., dust 
above sea spray), aerosol feature types that have been identified as marine in a given 
atmospheric column are not enough to carry out the analysis. Therefore, when de- 
termining the lidar ratio of SSA using Eq. (4), the algorithm only retains the data in 
which sea spray is the only type of aerosol present in the entire atmospheric column. 
To further reduce the uncertainty, we constrain the analysis to single layer profiles and 
remove profiles in which marine aerosol layers are vertically stacked within an atmo- 
spheric column. Therefore, the vertically integrated particulate attenuated backscatter 
T p is replaced by l~ p . Similarly, the column lidar ratio S p is reduced to S p for the re- 
mainder of the paper. Note also that all quantities discussed are particulate quantities 
(denoted by subscript p) and therefore, molecular scattering is removed using gridded 
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molecular and ozone number density profile data from the Goddard Earth Observing 
System Model, version 5 (GEOS-5) analysis product available from the NASA Global 
Modeling and Assimilation Office (GMAO) (Winker et al., 2009). Operationally, partic- 
ulate scattering is determined to be where the ratio of the CALIOP 532 nm scatter- 
ing profile normalised by the GEOS-5 molecular scattering profile is greater than one 

(jp- >1^. Errors associated with l~ p are discussed in Sect. 4. 

In this study, CALIOP and SODA nighttime data are collected globally for a period of 
three years from December 2007 to February 2011. The global data are then binned 
into 2° x 5° latitude and longitude, respectively, grid cells. The data is binned accord- 
ing to season and ocean surface wind speed (taken from the Advanced Microwave 
Scanning Radiometer-Earth (AMSR-E) observing system). To identify distinct features 
associated with the variability in sea spray aerosol lidar ratio over different parts of 
the oceans, the selected data is examined in relationship with other variables such as 
season, spatial location, and wind speed. 

For the data analyses conducted in this study, the single layer columns that have 
been identified by CALIPSO as clean marine aerosol were isolated and binned when 
the following conditions were met: (i) the vertical feature mask found one layer in the en- 
tire column, (ii) the vertical feature mask ranked the layer as type: aerosol and subtype: 
clean marine, (iii) the layer top was < 2 km, (iv) the relative error in r p due to random 
noise in molecular backscatter was < 50%, (v) the collocated SODA 5 km layer was 
composed of at least 70% shot-to-shot data (therefore increasing the signal to noise 
ratio), (vi) the total number of retrievals per 2° x 5° grid cell ranked above the first quar- 
tile of the grid cell frequency distribution (i.e., grid cells with frequency data in the lower 
25% of the distribution have not been included in the data analysis and are masked 
in the subsequent figures), and (vii) only nighttime data was used. Such strict quality 
controls, although reducing the total number of data points, considerably increase the 
reliability of the selected episodes. It should be noted that even after all the quality con- 
trol and quality assurance tests have been conducted, a large number (over 260000) 
of data points allowed a robust statistical analysis to be conducted. 
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To further minimise the effect of outliers on the estimated value, statistical medians 
were calculated for each grid cell. Regions where grid cells fail criteria (vi) indicate 
low confidence in the reported median value due to a low number of observations and 
are removed. This also ensures that global means are not heavily biased by outliers 
5 resulting from grid cells with a low number of retrievals. This procedure allows the 
removal of the non-physical positive skewness of the distribution observed by Josset 
et al. (2012). Despite such rigorous quality control criteria, readers should be cautioned 
when interpreting data near coastlines as the CALIOP scene classification algorithm 
may mistakenly identify continental pollution as clean marine aerosol (Burton et al., 
10 2012; Oo and Holz, 201 1 ; Schuster et al., 2012) causing an overestimation in the lidar 
ratio inferred from Eq. (4). Further discussion of error analysis is given in Sect. 4 below. 


3 Results 

3.1 Global distribution of retrieved AOD and lidar ratio 

Active detectors like CALIOP require knowledge of the lidar ratio for retrieval of aerosol 
is optical properties. Incorrect estimates of the S p values for a given aerosol type can 
lead to significant errors in the retrievals of particulate extinction and AOD. Past stud- 
ies using collocated CALIOP and MODIS retrievals have shown that, over the marine 
regions, CALIOP underestimates the AOD values relative to MODIS (Oo and Holz, 
2011). As MODIS data over the ocean has been extensively evaluated with numerous 
20 field campaigns (e.g. Levy et al., 2005), it was suggested that the primary source of 
discrepancy between the two sensors was the low value of the SSA lidar ratio used by 
CALIOP (Oo and Holz, 2011). Figure 1 shows seasonally averaged maps of CALIPSO 
and SODA SSA median optical depth at 532 nm and the differences between SODA 
and CALIOP retrieved AODs. White regions on Fig. 1 represent grid cells that were 
25 rejected by the data selection algorithm and have been removed from the subsequent 
data analysis. Inspection of Fig. 1 reveals considerable spatial and temporal variations 
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in sea spray AOD. Although the largest values of AOD seem to occur over regions with 
higher surface wind speed (i.e., the northern and southern oceans), elevated AOD val- 
ues can also be seen over the regions downwind from dust and/or pollution sources 
such as the mid-latitude North Atlantic Ocean and the Bay of Bengal (BoB) and over 
5 the major oceanic gyres. The former is believed to be just a retrieval artifact. Large 
disagreements between SODA and CALIOP reported AODs for these regions suggest 
that some dust/pollution aerosols might have been misclassified by CALIOP as sea 
spray. Higher S p values for dust and pollution compared to SSA would produce a higher 
AOD retrieval in SODA compared to CALIOP. Elevated AOD values over the oceanic 
10 regions with lower surface wind speed, on the other hand, could point to changes 
in sea spray particle size distribution to smaller sizes. Sub-micron sea salt aerosols 
(with particle diameter, D p < 1 pm) are believed to be characterised with larger lidar 
ratios than super-micron ones (e.g., Masonis et al., 2003; Oo and Holz, 2011). In gen- 
eral, Fig. 1 shows positive differences between SODA and CALIOP retrieved seasonal 
is median AOD values. Recalling that CALIOP retrieved extinction is the product of the 
prescribed lidar ratio and the measured column integrated particulate backscatter, pos- 
itive differences between SODA and CALIOP median AODs at 532 nm over most of the 
oceans suggest underestimation of the SSA lidar ratio prescribed in the CALIOP clean 
marine aerosol model. Figure 2 shows that over most of the ocean surfaces, the calcu- 
20 lated lidar ratio is higher than the default (S p = 20 sr) used in the CALIOP clean marine 
aerosol model. Global means and standard deviations for AOD and lidar ratio are given 
in Table 1 . CALIOP retrievals in this study cannot be directly compared to MODIS since 
we only use nighttime data. Nevertheless, SODA retrievals of AOD have been shown 
to agree well with MODIS (Josset et al., 2008) suggesting that the corrected lidar ratios 
25 will bring CALIOP retrievals close to MODIS data. Figure 2 also reveals that the value 
of the lidar ratio calculated using Eq. (4) changes considerably over different parts of 
the remote oceans, pointing to the variability in sea spray aerosol optical properties. 
It has long been known that meteorological and/or environmental factors and ocean 
chemical/biological composition influence SSA production, entrainment, transport, and 
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removal processes (Lewis and Schwartz, 2004), that can ultimately affect S p . More- 
over, due to atmospheric transport of SSA, satellite retrieved AOD values may also be 
related to the upwind processes. Despite the complexity of the mechanisms controlling 
SSA mass concentration over the oceans, surface wind speed has always been con- 
sidered as the major parameter governing the production, chemical composition, and 
life cycle of SSA (Lewis and Schwartz, 2004). Therefore, in the next section we will 
investigate the effect of wind speed on calculated temporal variability of marine aerosol 
lidar ratio. 

3.2 Wind speed dependence 

Numerous investigators have examined the effect of sea surface wind speed and sea 
state on marine aerosol optical properties. There are two production mechanisms for 
sea spray particles: bursting of bubbles at the water surface, and mechanical tear- 
ing of water drops (spume) from wave crests (for surface wind speeds U :o > 9 ms -1 , 
(Anguelova et al., 1999). Ocean bubbles are generated by the entrainment of air due to 
wave action. As bubbles rise due to their buoyancy, they burst at the surface producing 
SSA (Blanchard and Woodcock, 1957). In this study we have selected seven different 
wind speed regimes (see Table 2). The lowest wind speed regime, 0 < L/ 10 < 4ms -1 , 
was chosen to represent aerosols not generated via wind driven processes over the 
ocean. In general, ocean waves break at wind speed values above ~ 4ms -1 (initiating 
the white cap formation and bursting of the entrained bubbles) (Lewis and Schwartz, 
2004). Therefore, it has been suggested that below this threshold value, there should 
be a weak relationship between SSA optical properties and the surface wind speed 
(Kiliyanpilakkil and Meskhidze, 2011; Lehahn et al., 2010). Moreover, for such a low 
wind speed regime, most of the aerosols classified as clean marine by CALIOP are 
either produced outside the swath and then blown into the satellite field of view, or 
like in cases near coastlines, mistakenly identified as SSA. The highest wind speed 
regime, with l/ 10 > 15ms -1 , typically contributes a small fraction of CALIOP retrievals 
(Kiliyanpilakkil and Meskhidze, 2011) and is largely concentrated over the southern 
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ocean and in the northern Atlantic where the highest wind speeds are observed (Ben- 
tamy et al., 2003). Although CALIOP retrieval counts for sea spray aerosol in each 
2° x 5° grid cell are also influenced by the presence of clouds, Fig. SI shows the global 
distributions of CALIOP retrieval frequencies for different wind speed regimes. Figure 
5 3 shows the scatterplots for SODA and CALIOP retrieved AOD values for the wind 
speed regimes of Table 2. As expected, Fig. 3 shows that increases in wind speed 
are typically associated with higher values of SSA optical depth (note the center of 
the scatter distribution shifts to higher AODs for larger wind speed values). However, 
as the majority of the SODA AODs exist above the 1 : 1 line, this figure also indicates 
io the underestimation of CALIOP retrieved SSA optical depth values. When averaged 
over the entire globe, CALIOP retrieved clean marine AOD is roughly 32 % lower com- 
pared to SODA. According to Fig. 3 the largest discrepancies between SODA and 
CALIOP retrievals are observed at lower wind speed values. One easy explanation 
for this is greater chance for CALIOP misclassification over the oceanic regions where 
is long-ranged continental aerosols can contribute a larger fraction of the MBL particles 
(e.g., Blot et al., 2013). Terrestrial particles (e.g. mineral dust, anthropogenic pollution) 
are typically characterised by the larger lidar ratio values, leading to underestimation 
of the CALIOP retrieved AODs. However, measurements also show that changes in 
surface wind speed values can cause a considerable shift in the SSA size distribution. 
20 For optically active sea spray particles, the residence time decreases considerably with 
increasing size. Thus the aerosol population is increasingly controlled by the smaller 
end of the particle size spectrum as wind speeds decrease over the ocean (Hoffman 
and Duce, 1974). Conversely, as wind speed increases, fine mode aerosol volume size 
distribution changes slightly (with mixed trends), while the coarse mode volume size 
25 distribution exhibits a large and positive response to the increase in wind speed (Lewis 
and Schwartz, 2004; Smirnov et al., 2003). Such variability in SSA volume size distri- 
bution is expected to have an effect on the aerosol lidar ratio. As sub-micron sea salt 
aerosols are characterised with much larger lidar ratios than super-micron ones (e.g., 
Masonis et al., 2003; Oo and Holz, 2011), shifting SSA size distribution spectra to 
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smaller particles will cause an increase in total aerosol lidar ratio. Therefore, for clean 
marine aerosols, AODs and lidar ratios are expected to have opposite dependences 
on wind speed: high wind speed regions are characteristic of high AODs and low lidar 
ratios while lower wind speeds favor higher lidar ratios and lower AODs (Smirnov et al., 
5 2003; Sayer et al., 2012). 

Figure 4 shows that on average, the calculated aerosol lidar ratio is inversely related 
to the surface wind speed. According to this figure, aerosols retrieved in the wind speed 
regime 0 < U w < 4 ms -1 depict the largest variability in the lidar ratio as indicated by 
the spread of the distribution. As discussed above, aerosols in this regime likely include 
io both sea spray particles produced upwind and advected into the satellite field of view 
(with S p ~ 20 to 30 sr), as well as dust/pollution particles (with S p ~ 40 to 70 sr (Omar 
et al., 2009) that may have been misclassified by CALIOP as sea spray. As shown 
in Table 2, SSA lidar ratio distribution in this regime is characterised by the largest 
standard deviation (o = 1 7.4 sr) indicating that for the lowest wind speed values, a wide 
is range of sea spray particle sizes can be present over the ocean. Since for the wind 
speed values less than 4 m s -1 , the sea spray production is minimal, such large spread 
could also indicate that under low wind conditions there is greater probability for natural 
continental and human-induced pollution aerosols be miss-classified by CALIOP as 
clean marine. 

20 For the higher wind speed values (4 < U w < 15 ms -1 ) lidar ratio generally decreases 
with the increase in the wind speed and approaches the lidar ratios prescribed by 
CALIOP retrieval algorithms (i.e. , 20 sr) at the highest wind speed regime. According 
to Fig. SI, the most common wind values in CALIOP SSA retrievals over the ocean 
are in the 8 < L/ 10 < 10 ms -1 regime (26% of all available data) followed by the 6 < 
25 U w < 8 ms -1 regime (23% of all available data). For the higher wind speed regimes 
(U - io > 6 ms -1 ), surface winds play a decisive role in the determination of the lidar ratio 
(indicated by the narrow standard deviation, see Table 2). This is an important result as 
the distributions shown on Fig. 4 may help providing additional criteria for clean marine 
lidar ratio selection, yielding improved retrieval of sea spray AOD from CALIOP. 
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Analysis of data indicates that a mean lidar ratio of 26 sr is the most probable value 
that occurs for the majority of CALIOP retrievals over the oceans. This new lidar ratio 
reduces discrepancy between CALIOP-prescribed and SODA-derived lidar ratios from 
about 30% to 4%. Although the mean S p = 26 sr value for SSA proposed in this study 
considerably improves relationship between SODA and CALIOP retrievals for a wide 
range of ocean surface wind speed values, the added advantage of the current analysis 
is the ability to correct individual CALIOP retrievals of sea spray AODs. The correction 
to the lidar ratio can be prescribed using Eq. (4) as: 

1 _ e - 2T s = fE£(i _ e - 2T c) (5) 

s p,c 

where the subscripts C and S stand for CALIOP default (S p c ) and suggested (S p s ) 
lidar ratios, respectively. In future studies we intend to examine limited number of pa- 
rameters (in addition of wind speed) to construct a look-up table with a wind speed 
dependent, spatiotemporal distribution of SSA lidar ratios for use in CALIOP clean ma- 
rine aerosol retrievals. 


4 Assessing the validity and sensitivity of l" p 

The method used to derive the lidar ratio in this study depends on two parameters: 
the CALIOP integrated attenuated particulate backscatter (T p ) and the SODA aerosol 
optical depth (r p ). Uncertainties in both i" p and r p retrievals are expected to propagate 
through the calculations of the particulate lidar ratio. Josset et al. (2008, 2010a) in- 
vestigate the domain of validity for r p through an extensive calibration procedure and 
find that the SODA product, for retrievals at wind speeds between 3 and 10 ms -1 is in 
very good agreement ( R > 0.89) with MODIS AOD and has calibration errors less than 
1 5 %. Calibration errors in r p are expected to be even lower for nighttime retrievals used 
in this study (Josset et al., 2008). On the other hand, average uncertainty for CALIOP 
T p retrievals has not been examined previously and will be determined below. 
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Since ocean is the source of SSA, clean marine aerosol layers typically extend to 
the ocean surface. This makes it more difficult to determine molecular and particulate 
backscatter components of the signal separately using satellite measurements alone. 
To assess the uncertainty in lidar ratio introduced for the surface connected layers (i.e., 
layers whose bottom bound is defined as the ocean surface), here we estimate the 
error in CALIOP retrieved l~ p values. The total attenuated backscatter signal measured 
by the lidar consists of molecular and particulate components: 

P a u = (P P + PJe- 2Tp -e- 27 ™ ( 6 ) 

with subscripts m and p representing molecular and particulate quantities, respectively. 
From the definition of l~ p it follows that: 

z 

r P = j P p (z)e~ 2T pdz (7) 

o 

where the integration is from the surface to the top of the layer. (3 p is the particulate 

backscatter and e~ 2Tp accounts for the attenuation of the lidar signal by the particles. 
Substituting Eq. (6) into Eq. (7) gives: 

z 

r p = j {P ati e 2Tm - /? m (z)e“ 2T p)dz (8) 

o 

The molecular component of the signal in Eq. (8) can be derived from the GMAO 
modeled temperature and pressure profiles (Bloom et al., 2005). However, to solve 
this equation and determine the particulate attenuated backscatter value, particulate 
column integrated extinction is required. To get r p the CALIOP algorithm is using a pre- 
scribed value of the lidar ratio, making Eq. (4) circularly dependent on the lidar ratio. 
Here we estimate the error in CALIOP retrieved r D associated with the prescribed lidar 
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ratio by substituting the r p value from SODA. If the error is large, that would imply that 
the uncertainty in CALIOP prescribed lidar ratio would introduce sizable corrections to 
r p , making Eq. (4) unsuitable for the estimation of SSA lidar ratio. 

The relative error in f p can be defined as: 


r p ,s - r p ,c (e“ 2Tp ' c - e~ ZTpS ) • /g P m (z)dz 


2r n 


Error = 


p,c 


p,c 


( 9 ) 


where T p s and T p c are columnar integrated attenuated backscatter values for SODA 
and CALIOP, respectively. From the theoretical basis documents for CALIOP level 1 al- 
gorithms, the molecular backscatter is estimated as P m = where height depen- 

dent T (z) and P(z) profiles from the surface (1 000 hPa) to top-of-atmosphere (0.1 hPa) 
pressure levels were obtained from the GMAO Modern-Era Retrospective analysis for 
Research and Applications dataset. The molecular lidar ratio, S m is defined as 8n/3 
and C s is a constant equal to 3.742 x 10 -6 KhPa -1 m _1 (Hostetler et al., 2006). When 
considering all of the parameters, our analysis shows that the average error in f" p is 
approximately 1.5%. Compared to the systematic uncertainty in the SODA product 
(< 15%), the uncertainty in r p is much lower indicating that, on average, errors in l~ p do 
not dominate S p retrievals. Since an average discrepancy between CALIOP-prescribed 
and SODA-derived lidar ratios (~ 30 %) is more than an order of magnitude higher than 
uncertainty in f" p , we conclude that the uncertainty in the CALIOP column integrated 
backscatter has a minor effect on the Eq. (4) calculated lidar ratio. 

Furthermore, because in our study we use feature integrated products for a sin- 
gle aerosol layer, it is important to evaluate the relationship between T p and aerosol 
layer thickness (A Z). Figure 5 shows the normalised column attenuated particulate 
backscatter r p as a function of layer depth. For uniformly distributed aerosols through- 
out the column, r p is likely to be proportional to A Z. The spread of T p /AZ ratio is 
indicative of different amounts of SSA present in the column. Two limits of very high 
and very low A Z values are of particular interest. For example, strong reduction of the 
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r p /AZ ratio at the higher A Z values would indicate that the lidar signal is strongly 
attenuated throughout the layer reaching a sensitivity limit. On the other hand, con- 
siderable increase of the ratio for the thin layers may indicate contamination of the 
backscattered signal by the surface reflectance. According to Fig. 5 for the vast ma- 
jority of the data, signal attenuation and surface reflectance do not seem to be major 
issues for the surface connected layers, suggesting that the quality control algorithm 
described in Sect. 2.4 was sufficient to remove the majority of erroneous measures of 
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A new method is applied here that shows that it is possible to infer lidar ratios of sea 
spray aerosol over the ocean using two independent sources: the AOD from Syner- 
gized Optical Depth of Aerosols (SODA) and the integrated attenuated backscatter 
from Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). The proposed equa- 
tion calculates particulate lidar ratio for individual CALIOP retrievals of single aerosol 
layer columns as a correction to achieve the best agreement between SODA and 
CALIOP retrievals. The new method allows calculating sea spray lidar ratio and as- 
sessing its spatiotemporal variability and dependence on ocean surface wind speed. 
Analyses were carried out using CALIOP level 2, 5 km aerosol layer and collocated 
SODA nighttime data from December 2007 to December 2009. During the data analy- 
sis over 260 000 data points passed various quality-control and quality-assurance tests 
to reduce errors associated with the clean marine aerosol retrievals. The calculated li- 
dar ratios have been analyzed over the global ocean covering a wide range of wind 
speed and AOD conditions. Data analysis shows that over most of the ocean surfaces, 
the calculated lidar ratio is higher than the default lidar ratio of 20 sr used in the CALIOP 
clean marine aerosol model. The calculated aerosol lidar ratios are inversely related 
to the surface wind speed. Increase in mean surface ocean wind speed from 0 to 
> 1 5 m s _1 reduces the mean lidar ratios for sea spray particles from ~ 32 sr to ~ 22 sr. 
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Such reduction was explained by the shift in aerosol volume size distribution with the 
wind speed; however, it was also emphasised that future studies should explore the role 
of meteorological and/or environmental factors and ocean chemical/biological composi- 
tion for sea spray aerosol intensive properties. Our data analysis showed that changes 
5 in wind speed also affect the probability density function for sea spray aerosol (SSA) 
lidar ratio distribution. The largest standard deviation calculated for the lowest wind 
speed regime suggested that under low wind conditions, a wide range of sea spray 
particle sizes can be present over the ocean and there is greater probability for natural- 
continental and human-induced pollution aerosols to be classified by CALIOP as clean 
10 marine. Overall, our data analysis shows that an average value of 26 sr for SSA lidar 
ratio provides the best agreement between the SODA product and CALIOP retrieved 
global mean sea spray aerosol optical depth values. However, our study also shows 
large spatiotemporal variability in SSA lidar ratios, suggesting that a single constant 
value of the lidar ratio is not suitable for a wide range of sea spray aerosols and can 
is lead to large uncertainties at different locations and seasons. 

We have estimated the error in CALIOP retrieved column integrated attenuated 
particulate backscatter. Calculations suggest that the average error in particulate 
backscatter is more than an order of magnitude lower compared to the actual value. 
Data analysis also showed no clear indication for either approaching a sensitivity limit 
20 (due to strong attenuation of the lidar signal throughout the layer) or the contamination 
of the backscattered signal by the surface reflectance. Based on the conducted error 
analysis we conclude that the strict quality control criteria developed in this study is 
adequate to remove the majority of erroneous retrievals. 

Finally, even though calculations here were carried out for SSA, the technique used 
25 in this study is broad and can be used to infer lidar ratios of different species of atmo- 
spheric aerosols (i.e., mineral dust, biomass burning, etc.) advecting over the ocean. 
Because our data analysis shows that it is possible to derive a correction to the CALIOP 
prescribed sea spray lidar ratio, future studies should also consider conducting case 
studies over different oceanic regions to examine the possible effects of meteorologi- 
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cal parameters and ocean physiochemical/biological composition on sea spray aerosol 
lidar ratio. Classification (in the form of a look-up table) of spatiotemporal distribution 
and wind speed dependence of a limited number of parameters mostly affecting SSA 
lidar ratios, may lead to improved retrievals of AOD values over the oceans. 
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Table 1 . Seasonal means and standard deviations (in parenthesis) for 2° x 5° grid cell medians. 
The subscripts p, S, and C appended to r stand for particulate, SODA, and CALIOP, respec- 
tively, where r is the AOD. 
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Table 2. Means and standard deviations (in parenthesis) for 2° x 5° grid cell medians for various g 

AMSR-E wind speed regimes. The subscripts S and C appended to r stand for SODA and ^ 

CALIOP, respectively, where r is the AOD. 
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Fig. 1. Seasonal median AOD values from CALIOP and SODA (columns 1 and 2) and the 
difference (SODA -CALIOP) plot (column 3) for December-February (row 1), March-May (row 
2), June-August (row 3), September-November (row 4) plotted on a 2° x 5° latitude longitude 
grid. “No Data” is shaded white and is defined as grid cells failing quality control algorithm (see 
text for details). 
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Lidar Ratio, S p ( A = 532nm) 



Fig. 2. Seasonal lidar ratio for 2° x 5° latitude longitude grid cells. Seasons are arranged as (a) 
December-February, (b) March-May, (c) June-August, (d) September-November. 
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Fig. 3. Scatter plot of SODA (t s ) to CALIOP (r c ) AOD for each wind speed regime. Each 
crosshair indicates a grid cell median. The black solid and dashed lines indicate the 1 : 1 and 
least squares estimator lines respectively. Slopes (m) and correlation coefficients (R) are also 
reported. 
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Fig. 4. Probability density function of clean marine aerosol lidar ratio for selected AMSR-E wind 
speed regimes. The mean (/u) of each distribution is also reported. 
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Fig. 5. The normalised integrated attenuated backscatter as a function of the layer depth. The 
solid line shows the 3rd order least squares fit to the data while the dotted lines show ± 1 cr; the 
hatched area shows the layer depth data frequency: cross hatch between the 25th and 75th 
percentiles and straight hatch between 5th and 95th percentile. 
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